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ABSTRACT 

We develop an algorithm for the reconstruction of the two-dimensional mass distribution of a gravitational 
lens from the observable distortion of background galaxies. From the measured reduced shear 7;/(l — k) the 
lens mapping is obtained, from which a mass distribution is derived. This is unlike other methods where the 
convergence k is directly obtained. We show that this method works best for sub-critical lenses, but can be 
applied to a critical lens away from the critical lines. For finite fields the usual mass-sheet degeneracy is shown to 
exist in this method as well. We show that the algorithm reproduces the mass distribution within acceptable limits 
when applied to simulated noisy data. 

Subject headings: Cosmology: gravitational lensing; Cosmology: dark matter 
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1. INTRODUCTION 

The mapping of the distribution of matter in extended regions 
around rich clusters of galaxies, from the systematic distor- 
tions of background galaxies due to gravitational lensing, has 
increasingly become feasible and popular over the last decade. 
This shear is measured from the quadrupole moments of the 
images of the background galaxies in terms of the local 'polar- 
ization' of an image compared to its assumed intrinsic circular 
form. In order to remove the effect of the source ellipticity dis- 
tribution, the parameters are averaged over a large number of 
source galaxies. The two-dimensional mass distribution of the 
lens is reconstructed from the shear map away from the critical 
lines ('weak reconstruction'). Following the pioneering work 
of Tyson et al. (1990), the idea was quantitatively developed 
by Kaiser and Squires (1993, also Kaiser 1995 and Squires and 
Kaiser 1996) and Seitz and Schneider (1995, 1996). Since then, 
several interesting variants on this theme have appeared in the 
literature, e.g. using maximum likelihood cluster reconstruction 
(Bartelmann et al. 1996, Seitz et al. 1998), or using methods 
based on the variational principle (Lombardi and Bertin 1999) 
or maximum entropy (Bridle et al. 1998). 

Since the shear data contain complete information about 
the mass distribution in two independent fields 7i/(l — k) and 
72/(1 — k), it is possible to obtain several algorithms to estimate 
the mass distribution from the measured shear. The observed 
shear is available only on a finite grid. Due to the unknown in- 
trinsic distribution of the ellipticities of the source galaxies, and 
the effect of the distortion of the PSF due to observing condi- 
tions (seeing, tracking etc.), the measured shear is noisy. It is 
therefore desirable to develop new algorithms in the hope that 
they might be able to deal with the noise better than other meth- 
ods. 

In this paper we develop an algorithm (LM: Lens Mapping 
algorithm) for the lens mass reconstruction from the measured 
reduced shear, g, = 7/(1 - k). The method involves two steps. 
First, the lens mapping is reconstructed from the reduced shear, 
which can be done uniquely with the assumption that the lens 
mapping goes to identity far away from the lens. As the second 
step, we show in §3 that for a sub-critical lens, the surface mass 
density can be reconstructed completely from the derived lens 
mapping. In reality, however, the measured shear is available 



only in a finite region of the lens plane. In this case the lens 
mapping cannot be uniquely obtained and hence the LM algo- 
rithm exhibits a mass-sheet degeneracy, which is also present 
in other methods of mass reconstruction. In §3.1, we charac- 
terize this degeneracy in the context of the LM algorithm. The 
performance of this method with discretely sampled, noisy data 
is dealt with in §4, where we demonstrate the various features 
of the LM algorithm, by reconstructing the mass distribution 
for an analytically given shear field, and separately showing the 
effects of discrete sampling and of uniform noise (due to mea- 
surement error) on the input data. 

2. WEAK LENSING AND THE MASS-SHEET DEGENERACY 

For most practical applications of gravitational lensing the 
lens can be considered to be thin. Under the small angle ap- 
proximation the lens equation is given by 

y = x-W(x), (1) 
where the source angular position is denoted as y, the im- 
age angular position as x and the relativistic potential tp sat- 
isfies the equation k(x) = S/E cr ; t = |V 2 V>, where E cr ; t = 
{c 1 /AnG)(d s / d\d\ s ). Here d s ,dj and d\ s are respectively the an- 
gular diameter distance to the source, to the lens and between 
the lens and the source. The coordinates for source and image 
are small compared to unity and we can treat their components 
as Cartesian. It is convenient to define two planes: that contain- 
ing the source is called the source plane and that containing the 
cluster, the lens plane (or image plane). 

The coordinate differentials in the source plane and the cor- 
responding differentials in the lens plane are related through 
dyi = Mjjdxj and the inverse of the magnification matrix M is 
given by 



AT 1 =(!-«) 



-82 



(2) 



where g = g\+ig 2 = (71 + '72)/(l - «), 71 = 5(^,11 -^,22) and 
72 = "0.12 and the subscripts denote differentiation with respect 
to the two components of the image coordinates. From the 
quadrupole moments of the surface brightness of the images 
of the background galaxies the reduced shear can be measured 
unambiguously in the regions where k < 1 (for observational 
details see Kaiser 1999 or Bartelmann & Schneider 2000). 
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FIG. 1. — [a] The coordinate system used in this paper. We show the lens-plane images of two straight lines passing through (yi, V2) in the source plane; [b] An 
equally-spaced Cartesian grid (in the source plane) is lensed by a two-component mass model, both PICMDs with zero ellipticity and kg = 0.8 and 0.5 respectively 
(described in §3). The axes are marked in units of the core radius r c , assumed to be same for both the models. The separation between the centers of the two 
distributions is 10r c . 



The continuity of y(xi 7 x 2 ) implies dj k yi = d\-yu which, along 
with, eq. (2) gives 

K .' ( ';,' ~ K - A/', = G aj ~ > ( 3 ) 

where K = ln(l - n) and G" 1 = M~ l /(l - n). On multiplying 
(3) by the inverse of G" 1 , and taking the trace of the resulting 
equation gives us 

2 2 

V,ln(l-«) = EE [G'V'Aj-G^)] . (4) 
i=\ j=i 

This equation was first derived by Kaiser (1995). It is clear that 
replacing 1 - n on the left hand side with A(l - k), where A is a 
constant, does not affect the equation. Therefore any particular 
solution of this equation can be used to obtain a one-parameter 
degenerate family of functions, all of which satisfy (4). This is 
known as the mass-sheet degeneracy. 

3. THE LENS MAPPING ALGORITHM 

In general the lens equation (1) is a many-to-one mapping, 
and consequently the inverse of the lens equation x = x(y) has 
several branches, with no single branch completely specifying 
the lens mapping. However, if the lens is sub-critical, then the 
lens mapping is one-to-one everywhere. In this case both the 
mapping and its inverse are uniquely defined and completely 
specify the mass distribution. Here we show that in this case 
it is possible to obtain the lens mapping completely from the 
reduced shear, with the additional assumption that it goes to 
identity sufficiently far away from the lens. 

For a non-critical lens, open curves in the source plane are 
mapped to open curves in the lens plane. In particular, any in- 
finite straight line is mapped to an infinite open curve in the 
lens plane. In Fig. la, we schematically show the images of the 
two perpendicular straight lines passing through the source at 
(y\,yi)- The point of intersection of the two curves gives the 
position of the image (X\, X 2 )- 

From (2) we can obtain the equation which the image of an 
arbitrary infinite straight line in the source plane satisfies in the 
lens plane. In particular, we obtain the equations which map 



the coordinate grid lines of a Cartesian coordinate system in 

the source plane to the lens plane 

, n dxi g 2 , „ dx 2 g 2 ,,, 

dyi = i-> — = ; dy 2 = i-> — = — — . (5) 

dx 2 1-gi dx\ 

These are first-order ordinary differential equations and can 
be uniquely integrated through any point in the lens plane. If 
we consider the source position as the intersection of the lines 
y\ = constant and y 2 = constant in the source plane, then the im- 
age of this point will be the intersection of the images of these 
two curves in the lens plane. In this manner we obtain a one-to- 
one mapping from the lens to the source plane. Away from the 
lens singularities we can integrate the equations (5) to obtain 
the numerical solutions 

xi=X(x 2 ;x); x 2 = Y(xux) (6) 
where X (x 2 ; x) is the mapped curve for dyi = 0, passing through 
the point x m the lens plane, and Y (x\ ; x) is the mapped curve 
for dy 2 = 0, passing through the point x m the lens plane 
(see Fig. la). The auxiliary variable x m the above equations 
also explicitly represents the family of curves which the equa- 
tions (5) generate when integrated with the initial conditions 
X(xi) = Xi an d Y(x\) = X2. respectively, for the two equations. 

In Fig. lb we show the image of a Cartesian grid that is 
equally spaced in the source plane, lensed by a two-component 
model, where both clusters are represented by Pseudo Isother- 
mal Circular Mass Distributions (PICMD) with the same core 
radius r c , separated by 10r e along the x\ axis. The central val- 
ues of the dimensionless surface density for the two compo- 
nents are ko = 0.8 and ko = 0.5 respectively. Fig. lb illustrates 
the integral curves of (5). 

To obtain the source position corresponding to the image po- 
sition x we use the assumption that far away from the lens the 
lens mapping is y,- = x, (see Fig. la). Therefore the integrated 
curves should go to the original unperturbed source grid lines 
at large |x|. This is ensured by the following asymptotic condi- 
tions 

km X(x 2 ;x) = yr, hm Y(xux)=y2- (7) 

X2 — >±00 X\ — >±0O 

The lens equation can now be formally written as 

yi=X(oo;x); y 2 = Y(w;x), (8) 
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FIG. 2. — The reconstructed mass distribution, using the Lens Mapping Al- 
gorithm, for the model used in Fig. 1(b), shown as contours of equal surface 
density k. The axes are marked in units of the core radius r c . 
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FIG. 3. — The reconstructed mass distribution for the kq = 2 case. The re- 
construction works well outside the critical region. The "reconstructed" mass 
inside the critical lines shows up as noise, since the method fails there. 



where the auxiliary variables x are now interpreted as the im- 
age coordinates x. The trace of the magnification matrix which 
is given in terms of k as 

dX(oo;x) dY(oo;x) 

\ — + " =2[1- K (x)] , (9) 

OXl 0X2 

provides an estimate of k at any point. Fig. 2 shows the lens 
mass distribution thus reconstructed from the shear represented 
by Fig. lb, assuming that the shear is noise-free and is known 
everywhere. 

It might appear surprising that we have been able to obtain 
the lens mapping solely from shear. To understand how, let us 
recall that our assumption that the lens mapping becomes iden- 
tity far away from the lens removes the mass-sheet degeneracy 
We also note that it is essentially equivalent to the assumption in 
Kaiser's method that the mass density vanishes far away from 
the center, though quantitatively the notion of "far" is differ- 
ent in the two methods. For a finite field we will see that the 
degeneracy reappears. 

The above discussion applies to points away from the critical 
lines of the lens mapping. At these singularities, the eqs. (5) 
become non-integrable. This becomes apparent by writing the 
eqs. (5) in terms of 7,-s and k, 

dxi = 72 dx 2 = 72 

dx2 1 — k— 71 ' dx\ l-K + 71 
For a sufficiently smooth mass distribution, the potential and 
consequently 72 is non-singular everywhere. Since the Jaco- 
bian of the lens mapping, (1 -k) 2 --/ 2 , is positive far away from 
the lens, we conclude that (1 -k) 2 -^ 2 > "f 2 . From this inequal- 
ity it is clear that (1 -k) 2 -^ 2 can become zero only when the 
Jacobian ( 1 - n) 2 - 7 2 also becomes zero. Since ( 1 - k) 2 - 7 2 = 
implies that either 1 - k - 71 = or 1 - k + 71 = 0, we conclude 
that both the equations (5) cannot be integrated simultaneously, 
which implies that a correspondence between the lens plane and 
source plane cannot be obtained. Thus the LM algorithm can- 
not be applied in regions very close to the critical lines of the 
Lens mapping. 

3.1. Finite Field, Critical Lens 



The algorithm as described in the previous section is applica- 
ble to cases where the reduced shear is known everywhere. In 
practice, this information is available only in a limited region of 
the lens plane. We show here that the LM algorithm can be gen- 
eralized to this case. For a finite field, eqs. (7) can be applied 
only up to the edge of region in which data is available, and the 
correct limit cannot be evaluated. It is easy to see that due to 
the unknown magnification factor at the edge, these equations 
provide an incorrect measure of the derivatives in eq. (9). 

Let us consider the differentials of X and Y (as defined in 
eq. 6) at the edge x e , 

AX AT 
AX(x e 2 , X ) = — Ayi; AY(x e uX ) = -r- Ay 2 . (11) 
ZYyi Ay 2 

Since at the edge of the image the derivatives of X and Y are 
> 1 we obtain slightly higher values of Ay,. Substituting the 
value of the derivatives in eq. (11) we obtain 

AX(xl lX ) = [l-«(x e )][l- gl (x e )] Ay! , (12) 

AY (x\ , X ) = [1 - « (x e )] [l+ gl (x e )] Ay 2 • (13) 

This gives us a way of obtaining the exact value of the deriva- 
tives of the lens mapping, and thus completely removing the 
mass-sheet degeneracy. However, in reality we only know the 
reduced shear, and consequently the correction factors at the 
edges, only up to a factor of 1 — k. For finite fields we can 
impose the boundary condition k = to obtain the mass distri- 
bution. This corrects the generalized degeneracy up to a factor 
of(l-K). 

This algorithm can be applied to a critical lens as well. Since 
the algorithm depends on our ability to integrate equations (5) 
without hitting a singularity, we see that the method can be used 
for the region which lie outside the critical lines. In Fig. 3, 
where the model has kq = 2, we see that our algorithm works 
well outside the critical region. 

4. NOISY AND DISCRETELY SAMPLED DATA 

We now present a few examples to illustrate the application 
of the LM algorithm to noisy data. Since we wish to compare 
the original mass distribution with that reconstructed by the al- 
gorithm, we correct for the mass-sheet degeneracy in all the 
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FIG. 4. — The effect of the finite sampling of shear: the percentage error 
in the reconstructed mass distribution (using the LM algorithm), in the case 
where the reduced shear is known on a grid of size 20" X 20", and there is no 
noise. The mass is reconstructed using the LM algorithm on a 20 X 20 grid. 
The excess noise near the center of the distribution is due to the fact that the 
reduced shear is most rapidly varying close to the center, where the interpo- 
lated shear values are consequently noisier. Note that the rms noise outside the 
central region is less than 1%. 



FIG. 5. — The effect of the finite sampling and noise: the percentage error 
in the reconstructed mass distribution in the same case as the last one, but with 
noise with a constant signal-to-noise ratio of 10 added to the shear field before 
the application of interpolation. 



examples by the prescription given in § 3.1. We consider a sin- 
gle component lens modeled as PICMD. We simulate the signal 
by analytically calculating the reduced shear from a model. We 
illustrate the effect of finite sampling of data and the presence 
of noise with separate examples. 

We first consider the case where the reduced shear is known 
at all the points from an analytical model. We choose a lens 
with the central, dimensionless surface mass density given by 
kq = 0.7, and reconstruct the mass distribution by applying the 
LM algorithm on a 20 x 20 grid. The noise in the reconstruc- 
tion, which is very small everywhere (< 1%), can be entirely 
accounted for as arising from the inaccuracies in the evaluation 
of various integrals and derivatives required by the algorithm. 

To illustrate the effect of finite sampling of shear in the im- 
age plane, we consider a model with a lens velocity dispersion 
of cr v = 1100 km/s and with a core radius R c = 50 kpc. The 
source galaxies are assumed to be at the redshift z s = 1 and the 
lens is at zi = 0.2. We sample the reduced shear in pixels of 
size 20" x 20". Using a bicubic spline interpolation routine to 
evaluate the reduced shear at all the points, the mass distribu- 
tion is reconstructed using the LM algorithm on a 20 x 20 grid 
as before. The resulting percentage error in the reconstructed 
mass distribution (given by 100 x (k kc - k^)/ n tme ) is shown 
in Fig. 4. We note that the excess noise near the center of the 
distribution is mainly due to the fact that the reduced shear is 
most rapidly varying close to the center and therefore the inter- 
polated shear values are noisier at those points. This noise will 
vanish if the grid on which the shear is evaluated were made 
finer, but then in the real world, one is limited by the number 



of galaxies over which the value of the reduced shear has to be 
averaged, and we cannot use a substantially finer grid without 
increasing shot noise. 

As our last example we consider the same case as the last 
one, but add a uniform noise with a constant signal-to-noise 
ratio of 10 to the given shear field, before the application of 
interpolation. The noise in the reconstructed mass distribution 
(Fig. 5) can be seen to be around twice that of the noise in the 
shear field. 

5. CONCLUSIONS 

We have described a simple algorithm for mass reconstruc- 
tion of a cluster lens by directly obtaining the lens mapping 
from the measured reduced shear and taking its derivatives to 
compute the surface mass density n. 

The method works best for sub-critical lenses, where it can 
give the mass distribution at all points, but it can work for criti- 
cal lenses as well in limited regions of the lens plane, away from 
the critical lines. The algorithm is shown to have a mass-sheet 
degeneracy of the same type as exists in the other methods of 
reconstruction if the field in which the reduced shear is avail- 
able is finite. We have tested the algorithm on discretely sam- 
pled noisy (simulated) data and have found that it reproduces 
the mass distribution within acceptable limits. 
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